Scanning Tunneling Microscopy of a Luttinger Liquid 
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Explicit predictions for Scanning Tunneling Microscopy (STM) on interacting one-dimensional elec- 
tron systems are made using the Luttinger liquid formalism. The STM current changes with distance 
from an impurity or boundary in a characteristic way, which reveals the spin-charge separation and 
the interaction strength in the system. The current exhibits Friedel-like oscillations, but also carries 
additional modulated behavior as a function of voltage and distance, which shows the spin-charge 
separation in real space. Moreover, very close to the boundary the current is strongly reduced, 
which is an indication of the interaction strength in the system. 
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In the last two decades the interest in quasi one- 
dimensional physics has been spurred by experimental 
progress in constructing smaller and more refined struc- 
tures such as carbon nanotubes, atomic point contacts, 
and mesoscopic quantum wires produced by etching H or 
cleaved edge overgrowth [||. More recently it was even 
possible to produce single atomic chains by depositing 
gold on a silicon surface [|| . 

The theoretical foundation for describing interacting 
one-dimensional electrons was laid in the early 80's with 
the concept of a Luttinger Liquid (LL) Interest- 
ingly, the electron-electron interactions cannot ever be 
neglected in one dimension which makes those systems 
fundamentally non-Fermi-liquid like. The elementary ex- 
citations are described by separate spin and charge quasi- 
particles which move at different velocities || . The corre- 
lation functions are power-laws with exponents that are 
related to a single interaction constant g. 

Angle resolved photo-emission experiments have made 
some progress in identifying a possible signature of spin- 
charge separation in quasi one-dimensional compounds 
||,[3| . On the other hand, in mesoscopic wires most exper- 
iments focus on conductivity measurements jljfl]. How- 
ever, from those experiments it is very difficult to extract 
information about the fundamental interactions within 
the wire. To test the important theoretical concept 
of spin-charge separation in mesoscopic systems other 
methods must be considered. 

One difficulty in producing a good wire is the fact that 
even small impurity perturbations effectively cut the wire 
at low temperatures 0. However, such boundaries give 
rise to other interesting effects, which can even reveal the 
LL behavior as we will show below. One well-known im- 
purity effect in metals is the induced charge density fluc- 
tuation at twice Fermi wave-vector, the so-called Friedel 
oscillation. In the case of carbon nanotubes the Friedel 
oscillations have already been used to show the more 
complicated electronic structure |8|], which stems from 
a rolled up two dimensional graphite sheet. 

We now make predictions for an STM experiment 



along a simple quantum wire described by the LL model 
with an open end instead of leads i.e. electrons which have 
been confined to move in one dimension by clever gating 
or an appropriate deposit on a surface. We show that the 
spatial structure of the tunneling current reveals both the 
spin-charge separation and the interaction strength in the 
LL system. 

The STM current / is directly related to how many 
electron states are locally available in the LL system and 
in the tunneling tip. In particular, at position r and for 
a given tunneling voltage V we can write 



7(V» oc / duj N{uj,r)f(uj - V), 



(1) 



where N(ui,r) is the local density of states (DOS) in the 
LL with u measured relative to the Fermi energy, and 
/ is the DOS in the tip. We do not know the detailed 
properties of the tip, but we can assume that /(c<->) is 
smooth compared to the more singular structure of N(lo), 
so that / = const is a valid approximation. The DOS of 
the system is given in terms of the time-time correlation 
functions in the LL at position r 
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N(u, r) = — e^({^ a (r, t), (r, 0)})dt, 

^ 7r J — oo 
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where uj is measured relative to the Fermi energy. It is 
already well understood how to calculate the DOS as a 
function of energy u> with the LL formalism |^|| , but as 
described in this Letter it is the spatial structure as a 
function of distance r from a boundary that also explic- 
itly shows the spin-charge separation. 

The LL Hamiltonian describes one-dimensional elec- 
trons with short range interactions in the low- 
temperature limit below some energy cut-off A. In that 
limit, the dispersion is approximately linear and the 
electron-electron scattering rate can be taken as momen- 
tum independent. The high energy cut-off A is large com- 
pared to the temperature, but about one magnitude less 
than the bandwidth for typical lattice models. The elec- 
tron field ^^(x) is expressed in terms of left- and right- 
moving Fermions at the Fermi points ±Icf 
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We also define the Fermion currents J1/r ='■ ^l/r^l/r'-- 
Umklapp scattering is suppressed away from half filling, 
and remarkably all forward scattering processes can be 
described by expressing the spin and charge currents in 
terms of separate bosonic variables 4> c and <p s and their 
conjugate momenta IT C and Tl s 
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(n c/s + d x cj> c / s ) 

The LL Hamiltonian density is then written as 



(4) 



H 



[g-Hd x cb c ) 2 +gUl} + f[( dx cf> s ) 2 + ILt] (5) 



which describes two independent bosonic excitations for 
spin and charge separately. Due to SU(2) invariancc 
the spin boson is a free field, but the charge boson gets 
rescaled by the LL parameter g which is less than unity 
for repulsive interactions. 

After the linearization around the Fermi points in 
Eq. (||) we can write (omitting the spin indices a) 

(*(r,t)*t( r ,o)> = (V L (r,*)4(r,0)> + <V>rM)4M)> 
+ e i2fc -'-(^(r,t)4(r,0)} + e- i2fe ^(^(r,t)4(r,0)). (6) 

The last two terms carry a rapid oscillation with twice 
the Fermi wave- vector 2k ft reminiscent of Friedel- 
oscillations, while the first two terms are slowly varying. 
The tunneling current in Eq. ([l]) therefore has a rapidly 
oscillating Friedel part I osc and a uniform part / un j with 
a spatial dependence that is smooth compared to e l2k " r 

I(V, r) = I nai (V, r) + cos(2fc F r + <f>)I osc (V, r). (7) 

In a translational invariant system the left- and right- 
movers are uncorrelated (iPl^r) = an d we cannot ob- 
serve any spatial structure. A generic impurity, however, 
scatters left- into right-movers and the resulting corre- 
lation functions depend on the distance r from the end. 
Such boundary correlation functions have first been cal- 
culated for the spin-channel Jl(]] and later also for the full 
electron field fll}-|l3|| . We consider an open boundary at 
the origin r = of a relatively long system so that N(u>) 
is continuous (the other end of the system is far away 
and can be neglected for now). In that case we find for 
the uniform terms 
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and for the Friedel terms 
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where the exponents are given in terms of the interaction 
parameter g 

a = (I + g + 2)/8, b = (± + g - 2)/8, c = (± - g)/8. 

The short distance cut-off a ~ v/A is small compared 
to all other length-scales in the system, and we take 
a —> + in all following calculations. The correspond- 
ing expressions for left- and right-movers exchanged can 
be obtained by taking r — > — r. 

Let us first consider the uniform part of the current 
7 un i as determined by the analytic structure of Eq. (|^). 
A change of variables t' =t/r in Eq. (Q) and to' — rco in 
Eq. (Q) shows that the uniform current I un i is a function 
of the scaling variable rV 



(9) 



(10) 



For non-interacting electrons (g = 1, v c = v s = v) we get 
a = 5 , b = c = corresponding to a single pole of order 
one at t = ia in Eq. (||). The integration in Eq. (|||) gives a 
constant from the residue in the upper half plane, i.e. the 
DOS is independent of r and u>, and Eq. (^) simply gives 
^uni oc V for non-interacting electrons. However, even for 
small interactions the single pole splits into three singu- 
larities at t = ia and t — ±2r/v c + ia in Eq. (||). Close 
to the boundary the behavior of the Fourier transform 
in Eq. (||) is then dictated by the large time behavior of 
Eq. (||) and we find 

/uni OC r (l/9-9)/4 F (l/9+l)/2 for r < ( n ) 

i.e. a characteristic depletion as r — > for repulsive in- 
teractions g < 1. On the other hand far away from the 
boundary r 3> v c /V the behavior is dominated by the 
most divergent singularity and we find 
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y{l/ g +g+2)/A for r>Vc / V , 



(12) 



which is largely independent of r. However, the inte- 
gration of the deformed contour in Eq. (|^) around the 
branch-cuts of the weaker singularities also contributes, 
multiplied by a corresponding slowly oscillating "residue- 
factor" e ± 2lrw /' u <:. This results in an additional slowly 
oscillating contribution with cos(2rV^/w c ) that drops off 
with r^ 1 / 9-9 " 8 )/ 8 . The depletion with the slow oscilla- 
tions towards a constant current is depicted in Fig. (|l|) 
for <7 = 3/4 from doing the integrals numerically. None of 
the singularities in Eq. (||) depend on the spin velocity v s 
and therefore the uniform current / un i will not show any 
signs of spin-charge separation. Nonetheless the weaker 
singularities still give the characteristic slowly oscillating 
structure due to interactions. 

The effect of the last two "Friedel" terms in Eq. @ 
on the other hand will reveal the spin-charge separation. 
The amplitude 7 osc of the rapidly oscillating Friedel cur- 
rent has the same scaling form as in Eq. (0) but with an 



2 



uni 
6.0 



5.5 



dV 4 h 

2 




co [vjr] 

2 4 6 



[v c /V] 



10 



FIG. 1. The uniform current / un i in arbitrary units as a 
function of r for g — 3/4. The inset shows dI un i/dV as a 
function of u at a position r. 



even richer behavior for F. Already for non-interacting 
electrons we find a single pole at t — 2r/v + ia in Eq. (^) 
which results in N(uj,r) oc cos(2fc^?r + 2ruj/v). The in- 
tegration in Eq. ([!]) gives a strong r-dependence of the 
amplitude I OBC oc sin(— )/r ("amplitude modulation") as 
shown in Fig. (^) for g = 1. With interactions we now 
find four different singularities at t = ia, t = ±2r /v c + ia 
and t = 2r/v s + ia in Eq. (^|). The long-time behavior of 
Eq. (j|) is the same as for the uniform terms in Eq. (||), 
so that we get the same universal depletion for I osc as 
in Eq. ( |Tl] ) very close to the boundary r < v s /V. For 
larger distances from the boundary r ^> v c /V, however, 
the drop-off of the Friedel current is determined by the 
leading singularity in Eq. <M) 



r -(l/g+g+2)/8y(l/g+g+2)/8 k < g < 1 
r -(l+ 9 )/2y(l/3- 9 )/4 g< l 



(13) 



More importantly, the Friedel amplitude / sc has an oscil- 
lating superstructure from the residue factor of each sin- 
gularity. The strong amplitude modulation with sin 
has already been demonstrated for the non-interacting 
case, but the ratio of the velocities v c / v s can now be much 
larger than one. Therefore, we observe two separate spin 
and charge amplitude modulations of / sc with rV/v s 
and with rV/v c respectively (which are still smooth com- 
pared to the overall oscillation of 2kpr). This behavior is 
demonstrated in Fig. (||) for g = 3/4 and v c /v s = 5. The 
physical interpretation is that we observe the superpo- 
sition of all electron wave- functions in the energy range 
from to V which exhibits the spin-charge separation 
due to the interference from the boundary. 

We have shown that the spatial electronic structure of 
an LL indeed shows the signatures of the spin-charge sep- 
aration, but it is important to critically analyze to what 
extent this could be observed in a realistic STM experi- 
ment. As an example we consider the monoatomic gold 
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FIG. 2. The Friedel amplitude of the tunneling cur- 
rent in arbitrary units as a function of r for g = 1 and 
g = 3/4, v c /v s — 5. The inset shows dI osc /dV at a posi- 
tion r with an arbitrarily chosen phase e l2k F r — 1. 



chains on a slanted silicon surface Si(lll)-Au(5xl), which 
showed the signature of LL behavior in photo-emission 
experiments || . The spin-charge separation seems to be 
present for all excitations over the entire bandwidth (ca. 1 
eV), but to observe the particular STM structures that 
we predict here, the voltage has to be below the cut-off 
A which is about 0.1 eV for this system. The most eas- 
ily observable feature in an STM experiment is probably 
the depletion of the tunneling current as a function of 
distance near the boundary. The range of this depletion 
is given by r w vf/V, where vf ~ 6 x 10 5 ^ for the 
gold chains, so that for V < lOOmeV the range is at least 
r > 40A. Already from the shape of this depletion an ap- 
proximate estimate of g can be made from Eq. (11). Sec- 
ondly, it is important to analyze the Friedel oscillations, 
which is a more difficult task. For the gold chains the 
Fermi vector is fc_F = Tr/2a with an interatomic spacing 
of a = 3.83A, so that the 2kpr oscillations are commensu- 
rate with the lattice. This makes the Friedel oscillations 
easier to detect, but the small amplitude modulations 
in Fig. (||) may not show up very clearly. The Friedel 
amplitude is modulated by at most 30% by the charge 
waves, while the spin modulations are even smaller. The 
spin modulations are actually weaker for stronger inter- 
actions (about 10% for g ~ 1, but only 1% for g < 1/2). 
Nonetheless, even if only a small hint of those superstruc- 
tures shows up in an STM image, it should be possible 
to track those modulations for different voltages. This 
will change their period systematically and may make it 
possible to identify them unambiguously. Interestingly, 
the total modulation is much stronger without any spin- 
charge separation (g=l), so if it cannot be observed at all 
it would also be an indication of interaction effects. For 
larger voltages V > A it is more instructive to look at the 
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total density of electrons given directly by Eqs. (|[]|) in 
the limit t — > 0, which recovers the established drop-off 
of the Friedel terms with r - {1+9 ^ 2 but without 

any spin-charge modulations. 

The measurement of dl/dV as a function of voltage in 
an STM experiment (spectral mode) will give additional 
information about the local DOS. Close to the boundary 
dl/dV has considerable structure as indicated in the in- 
sets of Figs. The additional contribution dI osc /dV 
from the Friedel terms in Fig. (^) is smaller and depends 
on the choice of phase e l2kFr , but shows the separate 
spin-charge effects very clearly. The gold chains men- 
tioned above showed a surprisingly strong depletion of 
the spectral weight at the Fermi surface as a function of 
energy in photo-emission experiments ||. It would be 
interesting to see if STM experiments on the very same 
samples also show such a strong depletion of the DOS 
with a characteristic power-law V 2b . 

One approximation we made is that the DOS in the 
tip f(u>) is relatively smooth, which is not too restric- 
tive since the observed structures in Figs. can be 
reproduced for almost all forms of f(uj) as long as there 
is a sharp upper limit at V in the integral of Eq. (|l|) . It 
is also possible to average over both signs of the voltage 
to eliminate the DOS of the tip somewhat, since the an- 
ticommutator in Eq. (|^) gives the same singularities for 
t -» -t, i.e. N(to) = N{-uj). 

The distance from the tip to the sample is also im- 
portant, since the tip may influence the system and play 
the role of an impurity itself. The optimal distance can 
be determined in the actual experiment by operating in 
spectral mode somewhere in the bulk of the wire. If the 
dl/dV curve changes qualitatively as the tip-sample dis- 
tance is decreased (except for an overall scale), then this 
would be a clear sign that the tip influences the sample. 
In particular, if the tip starts to act as an impurity the 
dl/dV curve should show a stronger depletion with the 
power-law V^ -1 )/ 2 

So far we have only considered a perfectly reflecting 
impurity, because it is expected that any generic pertur- 
bation renormalizes to the open boundary fixed point [Q . 
However, in an intermediate range around a "weaker" im- 
purity the Friedel effects show a non-universal behavior 
described in terms of a form factor [jl4|]. The range of 
this so-called boundary shrinks to zero for a perfectly 
reflecting barrier. Also interesting are "active" impuri- 
ties that have a net magnetic moment or carry an elec- 
tric charge. In the presence of interactions those impuri- 
ties may be overscreened, i.e. the nearest electrons over- 
compensate for the impurity charge or spin and in turn 
get screened by the next-nearest neighbors, etc. This fi- 
nally results in a screening cloud which is also a 2k ft 
effect, since the impurity Hamiltonian H lmp induces a 
non-zero expectation value in the Friedel terms of Eq. (^) , 
i.e. {i> LYR^irap) 7^ 0. The presence of both backscat- 



tered and induced Friedel terms has recently been demon- 
strated for a two-channel Kondo impurity |l5|] . 

Finally, we must also consider the effect of a second 
boundary in a finite system at r = L. The correlation 
functions in Eqs. (|8|j9j) are then described by powers of 
sine functions sin ■ We ex P ec t that the spa- 

tial structure from the interference of the standing waves 
gives a similar picture as in Figs. (|l|^|) close to the bound- 
aries as long as V 3> ttv/L. However, a more dramatic 
finite size effect is a discrete spectrum N(ui) due to the 
appearance of (^-functions in Eq. (|^). This results in 
Coulomb-blockade-like charging steps in I(V), which can 
also reveal the spin-charge separation and the interaction 
strength in the LL []l3|| . 

In conclusion we have shown that the tunneling current 
has decaying Friedel-like oscillations in a range around a 
boundary, but additionally the Friedel amplitude carries 
a characteristic periodic modulation in real space which 
reveals the separate spin and charge parts of the electron 
wave-functions. The period of those modulations is a 
function of the tunneling voltage, which is assumed to 
be small. We also find a characteristic depletion very 
close to the boundary of both the Friedel current and 
the uniform current, which is immediately related to the 
interaction constant in the wire. 
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